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Abstract 

In this paper we make a systematic study for the semileptonic decays B/B s — > 
(77, r/, G)(l + l~, W, vv) by employing the perturbative QCD (pQCD) factorization approach, we 
firstly calculate the form factors FQ^J r q (q 2 ) and FQ^tp s (q 2 ) in the pQCD approach with the in- 
clusion of the next-to-leading-order (NLO) contributions, evaluate the form factors Fq^^ 3 (q 2 ) 
in the pQCD approach at leading order, and then calculate the branching ratios of the considered 
decays in the ordinary n-n' mixing scheme and the n-n'-G mixing scheme respectively. From the 
numerical results and the phenomenological analysis, we found that (a) the NLO contribution to 
the form factors F Q +J?(0) and F^^ Va (0) can provide ~ 12% enhancement to the leading order 
ones; (b) the NLO pQCD predictions for Br(B~ — > rf''l~v\) with I = (e, u) agree well with 
the data; (c) the pQCD predictions for the branching ratios of other considered decay modes 
are in a wide range of 10~ 5 to 1CP 11 , and generally have a moderate dependence on the mixing 
scheme; and (d) in the n-n'-G mixing scheme, the pQCD predictions for Br(B~ — > Gl~i>i) with 
I = (e, u, r) are at the level of 10 -5 . Our predictions will be tested by LHCb and the forthcoming 
Super-B experiments. 
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I. INTRODUCTION 



Pseudoscalar mesons rj and r{ are rather different from other light pseudoscalar mesons 
7r or K, not only for the large mass of the rf meson, but also for their mixing and the 
possible gluonic components inside 77 and 7/, as well as the so-called Kv^> puzzle. The two- 
body charmless hadronic B meson decays involving 77 W as final states and several different 
mixing schemes of 77 and rj', for example, have been studied intensively by many authors 
[ll-H] in the standard model (SM) and the new physics models beyond the SMp], The 
semileptonic B/B s meson decays, furthermore, are also studied frequently for example in 
Refs. [9l-ll2l]. and play an important role in improving our understanding about the nature 
of 77 and rj'. 

At present, many precision measurements about the branching ratios and CP violating 
asymmetries of relevant B meson decays, such as B — > Krj^ decays and B s — > J/^r]^'\ 
become available (l3|. and have been interpreted successfully within the framework of the 
SM, for example, in Refs. 0, by employing the pQCD factorization approach. As for 
the B ds — > (r],r]',G)(l + l~ ,lv,vv) semileptonic decays considered in this paper, only two 
of them have been measured experimentally [la], the corresponding world averages are[l~3| 

BR(B + r)l + v t ) = (0.39 ± 0.08) x 10~ 4 , 

BR(B + ->• r)'l + v t ) = (0.23 ± 0.08) x 10" 4 . (1) 

For most semileptonic decay modes considered here, obviously, the experimental mea- 
surements are still missing at present, and we have to wait for LHC measurements or the 
forthcoming super-B experiments. 

On the theory side, the early theoretical predictions as given in Refs. @, [l0| were 
Br(B ± — » rjlv) ps 4 x 10~ 5 and Br(B ± — > rj'lv) k 2x 10~ 5 , basically agree with the 
data in ([1]). In Ref. [16], Chen and Geng also found that the decays of Bd tS — > t]'l + l^ are 
sensitive to the possible flavor-singlet effects, although such contribution is small in size. 
As is well-known, the semileptonic decay modes B( s \ — > (77, t]')(l + l~, Iv, vv) are playing 
an important role in testing the SM and in searching for the new physics (NP) beyond 
the SM, and can also be used to extract the form factors of F 0j+ T (q 2 ) for B/B s — > r),r)' 
transitions. 

Analogous to a previous paper [l7j], we here will make a systematic study for the 
semileptonic decays B / B s — )■ (rj, rj', G)(l + l~ , W, vv), where G stands for the possible pseu- 
doscalar glueball. Up to the leading Fock state and leading order in a s , there are two 
Feynman diagrams as shown in Fig.l for the B — > r) qiS transitions and three Feynman 
diagrams as shown in the Fig. 2 for the B — > G transition. 



In a recent paper 18], Li, Shen and Wang calculated the NLO corrections to the B — > n 
transition form factors at leading twist in the factorization theorem. In this paper, 
based on the assumption of the SU(3) flavor symmetry, we extend the NLO results for 



the B — > 7r form factors as presented in Ref. [18( to the cases of B^ — > rj q ^ s directly. We 
also calculate the B^ — > r] g transitions form factors by employing the pQCD factorization 
approach [191 ]. 

This paper is organized as follows. In Sec. II, we present the ordinary 77-77' mixing 
scheme and the rj-rj'-G mixing scheme, and collect the distribution amplitudes for the 
r] q ,r] s and r] g . In Sec. Ill, we calculate the corresponding transition form factors, and show 
the differential decay widths for the considered decays. The numerical results and relevant 
discussions are given in Sec. IV. And Sec.V contains the conclusions and a short summary. 
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FIG. 1. Feynman diagrams of Br s ^ — > rj qjS transition. The representing the weak vertex 
where the final lepton pairs come from, and P denote rj q or rj s . 






FIG. 2. Feynman diagrams of 
while G stands for the glueball. 



(b) (c) 

rj g transition. The (^) representing the weak vertex, 



II. MIXING SCHEMES AND WAVE FUNCTIONS 



In order to check the mixing scheme dependence of the pQCD predictions, we will use 
two different mixing schemes in our calculation and compare the corresponding numerical 
results. The first mixing scheme is the conventional Feldmann-Kroll-Stech (FKS) r\ — 



mixing scheme [20|, in which the physical r\ and rj' can be written as 



cos « 
sin< 



- sin i 

COS(/) 



(2) 



where the flavor state r) q = (uu + dd) / y/2 and r/ s = ss, while is the mixing angle 20j . The 



relation between the decay constants (/^, /*,, /£) and (f q , f s , ) can be found in Ref. (!]. 



The chiral enhancement wig and have been defined in Ref. |21[ by assuming the exact 



isospin symmetry m q 



"0 



rrid- The three input parameters f q , f s and 4> in the FKS 



20 



mixing scheme have been extracted from the data of the relevant exclusive processes 
f q = (1.07 ± 0.02)/;, f a = (1.34 ± 0.06)/; and = 39.3° ± 1.0° with f n = 0.13 GeV. 

In the second mixing scheme: i.e. the T]-T]'-G mixing scheme as denned in Ref. [22 
the physical states rj, r/' and G are related to the flavor state r) q , rj s , and unmixed glueball 
state rj g , through the rotation 





(3) 



with the mixing matrix U (8, <p, 0g) 



U(9, 



cos (j) + sin 6 sin 0i Aq — sin <p + sin 9 cos d% Aq — sin^sin^c 
sin (ft — cos 9 sin 9{ A G cos — cos 9 cos 9{ A G cos 9 sin 0^ 
— sin 9i sin 0g — cos 9i sin 0c cos 0g 



(4) 



with the angle 



9 + 9i, 9i = 54.7° is the ideal mixing angle, and the abbreviation 



A G = 1 — cos 0c, where we have assumed that the r] 8 does not mix with the glueball. 
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For the r\ q and ri x meson distribution amplitudes, their two-parton quark components 
are defined as [23j, 241, 



( Vq (p 2 )\q*(z)q b p(0)\0) 



<%(P2)|^W4(0)|0> 



■—^=8 ab 
2^ Jo 



tee***" {hMthM*) + teWmg^s) 



+mg[ 75 (#-l)]^(a;)} 



+m^[ 75 (#-l)] /37 0' s (x)} , 



(5) 



where p 2 is the momentum of rj q , rj s or t/ 9 which will be given in next section, n = (1, 0, 0j_), 
and v = (0,1, 0j_) are the unit vectors in the light-cone coordinator system, N c = 3 is the 
number of colors and the expressions of the chiral enhancement scales mg and m s Q are [2l[ 



m q = 



m r 



m I 

11 r 2 2 i i 2 ■ _ 

— 1:1 = - — m„ cos <p + m„, sin — 
2m n 2m„ ' ' 



2 V^fs / > 



m; 



r/ , — mc) cos sin ( 



— — = [m 2 , cos 2 + m 2 sin 2 <\> ^— (m 2 / — m 2 ) cos < > sin < 

2m. 2m, ' ' ° r 



(6) 



In this work, we employ m q = 1.50 GeV and m s = 1.90 GeV [14j. And we adopt the 
distribution amplitudes 0q S (a;) (twist-2), <^ g (x) and 0|, s (x) (twist-3) the same as defined 
in Refs.@,[2fl 



2V2iVc 

fg(s) 
2^/2N~ c 



6x(l -x) 1 + ^ {s) C 2 3/2 (2x- 1) 

1 + (sOr/3 - ^ (s) ) C 2 1/2 (2a; - 1) 
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(7) 



- 3 \ Vscos + -p 2 q[s) {l + 6a* w ) \ Cl /2 (2x - 1) 



:{1-2X) 



2^2N, 
x(l - 10x + 10x 2 ) 



1 

2' 



1 + 6 ( 5?73 - -r^s - ^Pg(s) - g/tys)^ 



ff(«) 



(9) 



with the mass ratios p q = 2m q /m qq and p s = 2m s /m ss , and the Gegenbauer polynomials, 
C 2 1/2 (t) = i(3t 2 -l), C 2 3/2 (t) = ^(5t 2 -l), Q 1/2 (t) = i(3-30t 2 + 35t 4 ).(10) 

where t = 2x — 1. The Gegenbauer moments appeared in 0^ s (x), <^ iS (aO and 0* (x), are 
the following 25 



0.115, a^ (s) 



-0.015, ^ 9(s) = 0.013, ul q{a) 



(11) 



For the glueball rj g , its leading-twist distribution amplitude is defined by [24j, |27|, |28|] 

(12) 



(^(P2)|^(z)^ ] (0)|0) = i " r/fV ^ 



96 



n p p2 
n ■ P2 Jo 



1 dxe lxp2 - z 0G(X) 



x(l — x) 



where f Vg = f s for transition B s — >■ r] g , while f Vq = \/2f q for B — > r\ g transition, and with 
the notation 

Al(z)A b u] (w) = {A«(z)A b u (w) - At(z)Al(w))/2 (13) 

and the function 2^ 

G (x) = x 2 (l - x) 2 Bf ) C 5 1 /2 (2x - 1) , C? /2 (t) = 5t . (14) 

with factor Pf , Pf for the transition of B — >■ r] q , B s — )■ r/ s respectively. We adopt Pf = 
B S 2 = B 2 = 4.6 ± 2.5 0, Hi in this work. 



III. FORM FACTORS AND SEMILEPTONIC DECAYS 

In the B (for the sake of simplicity, we use B here to denote both B and B s meson) 
meson rest frame, with the mg stands for the mass of B meson, the momentum of B 
meson and the final pseudoscalar meson can be written as 

Pl = ^|(l,l,0 T ), p 2 = ^|p(0,l,0 T ), (15) 

with the energy fraction p = 1 — q 2 /m 2 3 carried by the final meson (here q = Pi — p 2 )- The 
light spectator momenta k\ in the B meson and k 2 in the final meson are parameterized 
as 

fcl = ( X l"^|> °> fc lT)> k 2 = (0, X 2P^|, km). (16) 

For B — )■ P transition, the relevant form factors Fo,+ ('? 2 ) an d Pr(g 2 ) have been defined 
for example in Ref. [30j with the relation F (0) = F + (0). For the sake of convenience, one 
usually define the auxiliary form factors fi{q 2 ) and f 2 {q 2 ) as 

(P(p 2 ) |6(0)7^(0) |S(p!)) = /i(g 2 W + / 2 (g 2 )p 2M (17) 
in terms of which F + (q 2 ) and Po^ 2 ) read as 

^(a = ^[/i(a + / 2 (? 2 )], (18) 



'•"«/-) = i/iiP) i + / 2 2 ] + U(g 2 )[l - / 2 J . (19) 



The amplitude for the B — > P transition form factors can be factorized into the convolu- 
tion of the wave functions for the respective hadrons with the hard-scattering amplitude. 
The wave functions are nonperturbative and universal. The factorization formula for the 



B — > P form factors is written as I23l. 31 



:v< v / ,.. ,2». j>, dz + d 2 z T dy d 2 y T 



{P(p 2 ) 16(0)7^(0) \B(p x )) =giC F N c / dx 1 dx 2 d 2 k 1T d 2 k 2T 



(2tt)3 (2tt)3 (20) 



xe- lfc2 ^(P(p 2 )|<( ?/ )^(0)|0)e*' 2 (0|6 a (0)gK^|P( J 9 1 ))T7^ 



In Ref. 18j, the authors derived the &>r-dependent NLO hard kernel for the B — > n 
transition form factor. We here use their results directly for B — >■ 7r transition processes, 
and extend the expressions of B — > tc form factors to the ones for Bi s \ — > rj q ^ transitions, 
under the assumption of SU(3) flavor symmetry. Based on the studies of Ref. [l8J, the 
hard kernel function H at the NLO level can then be written as the form of 



H = H^\a s ) + H^Xal) = [1 + F(x u x 2 , p,p f , V , (i)} H®(a a 



(21) 



with the NLO factor F(xi, x 2 , p, p/, r], (i) 



F(xx,x 2 ,p,Pf,p,(i) 



O'.s 



(p f )C F ( 21 p 



Ait 



m, 



13 
~2 



— In — - I In — 2~ + — In — + — In (xix 2 ) 



m 



B 



16 



1 



4 



ni 



+ - ln^ X\ + — In x\ In x 2 + ( 2 In — j- H — In p 1 In x\ 



7 



1 



4 



_1 m| 
2 d 2 



3 In 



B 



Cr 



+ 2 + 



101 „ 219 

7T 2 H 

48 16 



(22) 



where Ci = 25m s [18j, pf is the factorization scale and set to the hard scales t\ or U as 
defined in the Appendix, the renormalization scale p is defined the same as in Ref. [181 ]. 

In this paper, we use the same kind of distribution amplitudes for both B and B s 
meson as in Refs. 32 



[x, b) = N B x 2 {l — x) 2 exp 



1 / xm B 

2 \ W 



u 2 b 2 



(23) 



where ui, and ujb s are the shape parameters for B or B s transition respectively, and we take 
Ub = 0.40 ±0.04 GeV and u>b s = 0.50 ±0.05 GeV for B/B s transitions. The normalization 
factors Nb m are related to the decay constants /b w through 



B, 



dx<f)B (s) (x, b = 0) - . 



(24) 



In the transverse configuration b-space and by including the Sudakov form factors and 
the threshold resummation effects, we obtain the form factors /i,2((Z 2 ) and Fxlq 2 ) for 



6 



B s —> (Vs,Vg) an( i B ~~ ^ (.VqiVg) transitions: 



fi s ^ Vs {q 2 ) — 167rCVm| J dx\dx 2 J b\db\b 2 db 2 (f) B {xi, &i) 

x{[r (<^(z 2 ) - 0'(z 2 )) • /i i(xi, x 2 , h, b 2 ) - r x 1 pm 2 B <f)° '(x 2 ) 

xh 2 (x 1 ,x 2 ,b 1 ,b 2 )] ■ a a (ti) exp [-S Bv (ti)] 
+ [X! ( P (p a (x 2 ) - 2r ^(x 2 )) + 4r xi^(x 2 )] 

x/ii(x 2 ,xi,6 2 ,6i) • a s (£ 2 ) exp [-5 , Br ,(t 2 )] |, 



(25) 



(g 2 ) = 16rrCFiTi 2 B J dx\dx 2 J b\db\b 2 db 2 (f) B {xi,bi) 



x { [ [ (a; 2 p + 1) a (x 2 ) + 2r ^ Q - ^ 0* (x 2 ) - x 2 <f{x 2 ) + 3^(x 2 ) 

x/ii(xi,x 2 , 61,62) - r £im| (1 + x 2 p) CT (x 2 ) • ^2(^1, 2:2, 61, 6 2 ) j 
xa s (ti)exp [-SB»j(ti)] 

+2r ^ + 1^ P Or 2 ) • ^1(^2,^1,62,61) • a a (t 2 ) exp [- S Br} (t 2 )]} , 

F^ Vs (q 2 ) = SnC F m 2 B J dx x dx 2 J 6id6i6 2 d6 2 (l + r P )(j) B (x 1: b 1 ) 
x| roxiwigf (x 2 ) • h 2 (x 1 ,x 2 ,b 1 ,b 2 ) 



(26) 



+ 



a (x 2 ) - r o x 2 p (a; 2 ) + r ( - + x 2 ) 0'(x 2 ) + r o CT (a; 2 ) 



x/ii(xi,x 2 ,6i,6 2 ) • a s (*i) exp [-^(ii)] 

+2r (j) p (x 2 ) ^1 + ^ • h 1 (x 2 ,x 1 ,b 2 ,b 1 ) ■ a s (t 2 ) exp [-S Br) (t 2 )} }, 



(27) 



h {q) ~ 7E 



J dx\dx 2 J bidbi(j>B(xi,bi] 



4> G {x 2 ) 

x 2 (l - x 2 ) 



x\xj- h b (x 1 ,x 2 ,b 1 ) ■ a s (t b )exp[-S BG ](t b ) 



i 



+ / 6 2 * 2 



xix 2 p • h a (xi,x 2 ,b 1 ,b 2 ) ■ a s (t a ) exp[-S BG ](t a ) 



+[xt + (1 - xi)(l - x 2 )p] • h c (x 1 ,x 2 ,b 1 ,b 2 ) ■ a s (t c )exp[-S BG ](t c ) 



(28) 



V) 



dx\dx 2 / bidbi4>B(xi,bi] 



x 



x 2 (l - x 2 )p 



x\ 

p 



(j> G (x 2 ) 
x 2 (l - x 2 ) 

h b (x u x 2 , 61) • a a (t b ) exp[-S BG ](t b ) 



-J b 2 db 2 
+ 



+ x 2 ) ■ /i a (xi,x 2 ,&i,& 2 ) ■ tts(^a)exp[-S's G ](t a ) 



2(1 - x 2 ) - Xl (2 - x 2 - -) 

P . 



h c (x 1 ,x 2 ,bi,b 2 ) ■ a s (t c )exp[-S BG ](t c 



(29) 



Km 2 B C 2 F f g 

y/E 

(j> G {.x 2 ) 
x 2 (l - x 2 ) 

- / &2<^2 



+ 



dx\dx 2 \ bidbUl + r P )(j) B (xi,bi 



x - 2a; 2 ) • ^(^l, ^2, &i) • a s (t b ) exp[-S BG ](tb) 

+ x 2 ) ■ h a (xi, x 2 , 61, 6 2 ) ■ a s (t a ) exp[-5 , BG ](t a ) 



2(xi +x 2 — 1 



7' 



— 



h c (x 1 ,x 2 ,bi,b 2 ) 



xa s (t c ) exp[-S BG ](t c 



(30) 



where CV = 4/3, r = m 9(s) /m B , r P = (m Vq /m B ,m G /m B ) for B — >■ (rj qi r] g ) transition, 
and rp = (rn v Jm Bg ,m G /m Bs ) for £> s — >■ (r) s ,r]g) transition. We employ m % = 0.18 ± 
0.08GeV, m Vs = 0.69GeV [Hj and m G = 1.4GeV Q in this work. For B r) q transition 
form factors, there is a factor of l/y/2 at the front of each expression of fi, 2 (q 2 ) and Fx(q 2 ), 
because there is only the u or d quark component of the r\ q involved in B — > rj q transition. 
And f g = \/2f q ( f g = f s ) for B — > rj g (B s — > r] g ) transition. The hard functions ^1,2,0,6,0 
the scales £i, 2 , a ,& jC and the Sudakov factors exp[— S Brj (t)] and exp[— S BG (t)} are given in 
the Appendix A. 

One should note that fi )2 (q 2 ) and -Ft(<? 2 ) as in Eqs. (I25H27P do not include the NLO 
corrections. In order to include the NLO corrections, the a s in Eqs. ( l25]|2~Tj) should be 
changed into a s -F(xi, x 2 , p, pj, p, Ci), with the NLO factor F(xi,x 2 , p, pj, p, Ci) as defined 
in Eq. (|2"2|) . As for Bt s \ — > i] g transitions, only the leading order expressions are available 
here. 

For the detail discussion of the effective Hamiltonian of the b — > ul~vi and b — > svv 
transitions, one can see Ref. 17|. The effective Hamiltonian of the b — > sl + l~ can be 
found in Appendix B. The differential decay widths of the b — > ul~i>i, b — > sl + l~ and 
b — >■ svv transitions are given as 17 



dT(b -»■ ulv x ) 
dq 2 



G 2 F \V ub \ 2 q 2 - m] 
1927r 3 m| 



(q 2 



m 



2\ 2 



ml 



ml 



— m 2 p 



(q 2 ) 2 y q 2 y 4g 2 
x { (m 2 + 2q 2 ) [q 2 - (m B - m P f] [q 2 - (m B + m P f] F 2 + (q 2 ) 

+3m 2 (m 2 B -m 2 P ) 2 F 2 (q 2 )}, (31) 
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dT(b^sl+l-) G 2 F a 2 em \V tb \ 2 \V*\ 2 ^\(^) q 2 -Am 2 1 



dq 2 



512m 3 B ii 5 



3q 2 



x 



6m 2 \C 10 \ 2 (m 2 B - m 2 P ) 2 F 2 {q 2 ) + \C w \ 2 (q 2 - Am 2 ) \{q 2 ) F 2 + (q 2 ) 



+ (g 2 + 2m 2 )A(g 2 



C eff F+{q 2 ) + 2Cr ff {m h -m s )F T { l h 



tub + mp 



(32) 



dT(b — y svv) 
dq 2 



\G 



6->s|2 



\ 3 / 2 (m 2 B , nip, q 2 ) 
9Qm 3 B n 3 



2\|2 



(33) 



where Gp = 1.16637 x 10~ 5 GeV~ 2 is the Fermi-coupling constant, V u b,V t b and V ts are the 
CKM matrix elements, mi is the mass of the charged lepton, X(q 2 ) = (m 2 B + m 2 P — q 2 ) 2 — 
4m 2 B m 2 P is the phase-space factor and a em = 1/137 is the fine structure constant. The 
factor 3 in Eq.( l33l) arises form the summation over the three neutrino generations. The 
differential decay widths of b — >■ dl + l~ and b — > dvv can be obtained from the replacement 
of \V^\ ->• |V£|, m s ->• m d and \C h ^ s \ \C b L 

= \ct d \ x \v td /v*\. 



in 



" — > \C b T ^ d \. The expression of \C b £* s \ can be found 



171 ]. and we employ the relation \C L 



IV. NUMERCIAL RESULTS AND DISCUSSIONS 



In the numerical calculations we use the following input parameters (the masses and 
decay constants are all in unit of GeV) [lH, 20], 



A 



g" 4) = 0.287, f B = 0.21, f Bs = 0.23, 
ms± = 5.2792, mpo = 5.2795, m B o = 5.3663, 
t b o = 1.525 ps, r B o = 1.472 ps, m T = 1.777, 



Tp± = 1.638 ps, 
m b = 4.8, 



m w = 80.4, m t = 172, 



0.548, m v , = 0.958. 



(34) 



For relevant CKM matrix elements we use |V t6 | = 0.999, \V ts \ = 0.0403, \V td /V t 



0.211 13 



ts 



A. Form factors in the pQCD factorization approach 



As we have explained in [171 ]. the pQCD predictions for the considered form factors 
Fo(q 2 ), F + (q 2 ) and _FV(g 2 ) are reliable only for small or moderate values of q 2 . For the 
form factors in the larger q 2 region, one has to make an extrapolation for them from the 
lower q 2 region to larger q 2 region. For the form factors of B/B s — > {r) q ,r) s ,n g ) transition, 
we make the extrapolation by using the pole model parametrization 



*l(0) 



a(q 2 /m 2 B ) + b(q 2 /m B ) 2 



(35) 



where Fi denotes any function among -Fo i+i t, anc ^ a ' ^ are ^ ne constants to be determined 
by the fitting procedure. 
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TABLE I. The pQCD predictions for form factors .Fo,+,t(0) and the parametrization constant 
a and b for B — >• r] q and B s — > r] s transitions at the LO and NLO level. 







F(0) 


a 


b 


^0 


LO 
NLO 


0.17 ±0.02 ±0.02 ±0.01 
0.19±g;g| ± 0.02 ± 0.01 


0.57 ±0.02 
^2 +0 - 04 

u - oz -0.03 


-0.02 ±0.01 
0.00 ±0.01 


r + 


LO 
NLO 


0.17 ±0.02 ±0.02 ±0.01 
0.191:^1 ± 0.02 ± 0.01 


1.47 ±0.02 
1.43 ± 0.03 


0.44 ± 0.02 
n41 +o.04 

u - 4i -0.03 


T 


LO 
NLO 


0.15 ±0.02 ±0.01 ±0.01 
0.17 ±0.02 ±0.02 ±0.01 


1.54 ± 0.02 
1.50 ± 0.04 


o ™+°- 03 

u - ou -0.02 

n 47+0-05 

u - 4 '-0.04 


^0 


LO 
NLO 


0.27±g;ga ± 0.04 ± 0.01 
0.31 ± 0.04 ± 0.01 


0.57 ±0.02 
0.53 ±0.03 


-0.01 ±0.01 
0.00 ±0.01 


r + 


LO 
NLO 


0.27t^3 ± °- 04 ± °- 01 
0.31 to^t ± 0.04 ± 0.01 


1.52 ±0.02 
1.48 ±0.03 


0.48 ± 0.02 
0.46 ± 0.03 


rpBs^-ris 
r T 


LO 
NLO 


0.27 ±0.04 ±0.04 ±0.01 
0.31±g;g| ± 0.04 ± 0.01 


1.57 ±0.02 
1.53 ±0.03 


0.53 ±0.03 
0.49 ± 0.04 



TABLE II. The pQCD predictions for form factors F 0i+i t(0) (in units of 10~ 2 ) and the 
parametrization constant a and b for B — )• ij g and B s — > r\ g transitions. The largest error 
of -Fo,+,t(0) comes from the uncertainty of B<i = 4.6 ± 2.5. 





F(Q) 


a 


b 


TpB^g 

^0 

r B^r, g 
r + 

r T 


0.14±g;g? ± 0.01 ± 0.08 
0.14±g;gf ± 0.01 ± 0.08 
0.10 ±0.01 ±0.01 ±0.06 


1.20 ±0.01 

1.85 ±0.00 
1 q«+o.oi 


u - 2 °-o.oi 
0.86 ±0.00 

1.03 ±0.03 


rpBs^T/g 
^0 

rpBs^r/g 
r + 

T 


0.11 ±0.01 ±0.01 ±0.06 
0.11 ±0.01 ±0.01 ±0.06 
0.08 ±0.01 ±0.01 ±0.04 


1.15 ±0.03 

i o 4 +0.01 

i qo+o-o 3 


0.24 ±0.03 

n OO + 0.02 

U - 8ci -0.03 

n Q4+0-05 

u - J ^-0.03 



In Table |T]we list the LO and NLO pQCD predictions for the form factors -F o (0), F + (0) 
and Ft(0), and the corresponding parametrization constant a and b for B — >■ r] q and 
B s — > i] s transitions. The first error of ^,+^(0) in table U comes from the uncertainty 
of ojb = 0.40 ± 0.04 GeV or wg s = 0.50 ± 0.05 GeV, the second error comes from the 
uncertainty of f B = 0.21 ± 0.02 GeV or f Bs = 0.23 ± 0.03 GeV, while the third error is 
induced by a^ 3 = 0.115 ± 0.115. The errors from the uncertainties of Vt s or |V^/Vj s | 
and f q = (1.07 ± 0.02)/,,- or f s = (1.34 ± 0.06)/,,- are very small and have been neglected. 
The error for parameters a and b in Table [J comes from the uncertainty of a^^" , others 
are negligibly small. 

In Table [TTJ we list the pQCD predictions (in units of 10~ 2 ) for the form factors 
-F (0),F + (0) and Fjn(0) and the corresponding parametrization constant a and b for 
B — >■ T] g and B s — > rj g transitions. The source of the errors is the same as B — > r] q 
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FIG. 3. The pQCD predictions for the g 2 -dependence of F 0j+j t(9 2 ) f° r B — ^ rj q and B s — > r/ s 
transitions, where the solid curve, the dots curve and dot-dashed curve shows total, LO and 
NLO contributions respectively, and the shaded band shows the total theoretical error of the 
pQCD predictions. 




01234567 01234567 01234567 

q 2 (GeV 2 ) q 2 (GeV 2 ) q 2 (GeV 2 ) 



B -> Tl 

s 'g 













0.4 

0.3 

\ 0.2 
u_ 

0.1 
0.0 



- B -> Tl 











1 2 3 4 5 6 7 
q 2 (GeV 2 ) 



1 2 3 4 5 6 7 
q 2 (GeV 2 ) 




FIG. 4. The pQCD predictions (in units of 1CT 2 ) for the c/ 2 -dependence of Fo i+j r(g 2 ) for 
B —7- rig and B s — > r] g transitions, where the shaded band shows the total theoretical error of 
the pQCD predictions. 



and B s — > t] s transition except the third one from the uncertainty of B 2 = 4.6 ± 2.5. The 
errors of a and b here are induced by B 2 = 4.6 ± 2.5. 

In Figs. ([31 HD we show the pQCD predictions for the g 2 -dependence of the B — » rj qi 
B s — )• r] s , B — >■ r/g and B s — > rj g transition form factors F (q 2 ), F + (q 2 ), -Pr(g 2 ) , where 
the solid curve, the dots curve and dot-dashed curve shows the total, LO and NLO 
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contributions respectively, while the shaded band shows the total theoretical error of the 
pQCD predictions. 

From the numerical results as listed in Table HJ Table [III and the q 2 - dependence of the 
i*b,+,T as illustrated in Figs. ([3j HJ), we find that: 

(i) Although F (0) = F + (0) by definition, they have a similar but a little different q 2 
dependence, as illustrated by Fig. El The NLO contribution to the form factors of 
B(s) Vq(s) can provide about 12% enhancement to the leading order predictions. 
Since the B/B s — > rj g transition form factors are proportional to B 2 , they have a 
large theoretical error due to the large uncertainty of B 2 = 4.6 ± 2.5, as illustrated 
by Figs]3] and H 

(ii) For B/B s — > rj g transitions, the corresponding form factors Fo i+i t(0) are at the 10~ 3 
level, smaller than those for B — >■ r] q and B s — > rj s transitions by about two orders 
of magnitude, i.e. 



F t B ~^"(0) F? s ^ a (0) 



0.01, (36) 



where i = (0, +, T). In Ref.[34j, the authors calculated the form factors of B — > G 
and B ->• S and found that F?^ G (0) / Ff^ s (0) ~ 0.1 instead of 0.01. The major 



reason for this difference is that the glueball "G" in Ref. [34j is a scalar instead of 
a pseudo-scalar in this paper. For a pseudo-scalar glueball, the leading Gegenbauer 
moment in its wave function is zero, only the higher order Gegenbauer moment are 
non zero and contributes, so the form factor of B meson to pseudo-scalar glueball 
transition is much smaller than the one for B to a scalar glueball. After all, the 
contributions from the Bj B s — >• r] g transitions are indeed negligibly small. 

(iii) For the B — > r] q , B s — > r] s and B/B s — > i] g transitions, the pQCD predictions for the 
form factors Fq :+j t(0) agree with the values in literature. For F^~^ riq (0), for example, 
the NLO pQCD prediction is F+^ V9 (0) = 0.19 ± 0.04, which agree with the value 



as given in Refs. [24j, |33| when we take the difference in definition into account. For 
F+~* V9 (0), for instance, the pQCD prediction is F^^_{0^ = (0.14 ± 0.08) x 10" 2 , 



which also agree with the value as given in Refs. [241 . |27 



B. The branching ratios of the considered decays 

By using the relevant formulaes and the input parameters as defined or given in previous 
sections, it is straightforward to calculate the branching ratios for all the considered 
decays. 

We firstly consider the FKS mixing scheme: i.e. we assume that there is no rj g com- 
ponent in 77W meson. In Figs.fJMS]), we show the differential decay rates dT/dq 2 for the 

1 If wc consider the constant difference l/4iV c = 1/12 in definition and the second factor of 2 from the 
two gluons crossed diagram absorbed by the rj g wave function, the form factor F^^ Vs (0) = 0.044 as 
given in Ref. Q will be transformed into F^ Vg (0) = 0.044/(47V c )/2 = 0.18 x 10~ 2 , which agrees well 
with the pQCD prediction given here. 
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FIG. 5. The g 2 -dependence of the differential decay rates dT/dq 2 for the decays B — > 
ri^e~v e ( y [i~Vy) and B~ — > t]^t~v t . 

decay modes corresponding to the B/B s — > ry-'> transitions. By making the numerical 
integrations over the whole range of q 2 , we get the numerical results for the branching 
ratios. For b — > u charged current processes the pQCD predictions for the decay rates are 
the following 

Br(B~ 7,1-*) = (0.41±^ 2 ± 0.081J8S) x x W~\ (37) 



V, 



|2 



Br(B- -+ W -v T ) = (0.24±g^±85| ± 0.02) x x 10" 4 , (3 



0.0038 



Br (B- = (0.20±°;°4 ± 0.04 ± 0.02) x '^ |2 x 10" 4 , (39) 



038| 2 

r |2 
ufe | 

0.0038I 2 



5r(jB - ^ r -i7 r ) = (0.10j£jg ± 0.02 ± 0.01) x ffi 1 . x 10 -4 j (40) 



where the first error comes from the uncertainties of Ub = 0.40 ± 0.04, the second error 
comes from the uncertainty of fs = 0.21 ±0.02, the third one is induced by the variations 



of a 2 = 0.115 ± 0.115, and we adopt \V u b\ = 0.0038 [17J in the calculation. For both 



B~ — > r,l~Vi and B~ — > r,l~Vi decay modes, it is easy to see that the pQCD predictions in 
Eqs. (I37|39l) agree well with the measured vales as given in Eq. ([I]). The pQCD predictions 
as given in Eqs. ( I38"f4*0"]) will be tested in LHCb and forthcoming super-B experiments. 

For other neutral current processes the pQCD predictions are all listed in Table IHII 
The first error comes from the uncertainties of cj& = 0.40 ± 0.04 or ub s = 0.50 ± 0.05, 
the second error comes from the uncertainty of fs = 0.21 ± 0.02 or fs a = 0.23 ± 0.03, 
the third one is induced by the variations of a^' r,s = 0.115 ± 0.115. From the numerical 
predictions in Table III, one can see that Br(B° s -> tjvv) » Br(B° s -> rfuv) » 1.7 x 10~ 6 , 
while the branching ratios of other decays are at the level 10~ 7 — 10~ 9 . It is in fact a very 
hard job to observe these very rare decays even if it is not impossible. 
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FIG. 6. The g 2 -dependence of the differential decay rates dT/dq 2 for the decays B° / 'B® — > 
Tf('> (e + e~ , fi + [i~ , t + t~ , vv). 

In the FKS mixing scheme, there is no mixing between the pseudoscalar glueball and 
the r) q and r] s . The branching ratios Br(B° — > Gl + l~) with / = (e,/i, r) will be smaller 
than 1CT 11 , and can not be measured even at future super-B experiments. 

In the second rj-rj'-G mixing scheme, the physical states rj and rj' along with the pseu- 
doscalar glueball G are related to the flavor states rj q , r] s and rj g . From Eqs.([3]|l]) it is easy 
to find the following relations 

\r]) = (cos (ft + sin 9 sin Q { A G ) \r] g ) — (sin <p — sin 6 cos 0j A G ) \r] s ) — sin 9 sin G \rj g ) , (41) 
|?/) = (sin0 — cos^sin^iAc) \r] q ) + (cos0 — cos^ cos^Ag) \rj s ) + cos 9 sin <pc |?7 3 ) ,(42) 
\G) = -sin6»iSin0 G \r] q ) -cos6»iSin0 G \q a ) +cos0 G \q g ), (43) 

Since the form factors of B/B s — > r\ g transitions are two orders smaller than those of 
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TABLE III. The pQCD predictions for the branching ratios of the considered decays with I 
(e, /i), the errors comes from the uncertainties of wg, the /b, > and a^" 71 " respectively. 



Decay modes 



Br(B° 
Br(B° 
Br(B° 



Br(B° 
Br(B° 
Br(B° 



Br{B° s 
Br(B G s 
Br(B° s 



Br(B° s 
Br(Bg 
Br(B° s 



n 



rjl 

T)T ' T 

rjuv) 



rfl+V 

Tj'lSu) 



r,l+l-\ 
rjisis) 



rfi+r 

rj'isu) 



pQCD predictions 



(0.48 
(0.98 
(3.84 



0.14+0.10+0.05 
0.10-0.09-0.04 
0.28+0.19+0.08 
0.20-0.18-0.06 
1.11+0.76+0.33 
-0.83-0.70-0.32 



) X 10- 8 

X 10~ 9 

x 10~ 8 



(0.24 
(0.25 
(1.83 



0.07+0.05 

0.05-0.04 

+0.07+0.05+0.02 

-0.05-0.04-0.01 
0.53+0.37+0.m v 1fl 
0.40-0.33-0.15^ x 1U 



±0.02) x 10~ f 
x 10~ 9 

8 



(2.07 
(0.45i 
f 1.62 



+0.65+0.57+0. 10\ v ln -7 
-0.51-0.50-0.09^ A ±u 

g;i^:i?±0.02) xio- 7 

+0.54+0.45+0. 11 \ v ln -6 
-0.38-0.39-0.10y x ±u 



(2.18; 
(0.27; 

d.7i 



-0.73H 
-0.52- 
-0.09 



0.61 
0.53 



±0.15 x 10 



-7 



-0.07 



±0.07 ±0.01 x 10" 



-0.57H 
-0.41- 



0.47 
0.42 



±0.12) X 10" 6 



B/B s — » r] q ,r] s transitions, the contribution of the third part at the right hand side of the 
Eqs.( l41~H43|) can be neglected safely, and then we have 

1 77) = (cos0 + sin 9 sin Q { A G ) \rj q ) — (sin0 — sin 6 cos 6i A G ) |r/ s ), (44) 
\r]') = (sin — cos 8 sin ^ A G ) \rj q ) + (cos <fi — cos 9 cos Q{ A G ) \rj s ) , (45) 
\G) = - sin 6i sin 0g|%) - cos^ sin0 G |r/ s ), (46) 

In the Table IIVI we list the pQCD predictions for the branching ratios in the rj-rj'-G 



mixing scheme with the choice of = 43.7°, <p G = 33° 14], or = 40°, G = 22° [22 
From Table IIVI we observe that 

(i) The pQCD predictions for the branching ratios in the r]-T]'-G mixing scheme are 
close to but generally smaller than those of the corresponding decay modes in the 
FKS mixing scheme, except for the decay modes of B° s — > r// + / _ (r + r _ , vv). 

(ii) The pQCD predictions with the choice of = 40°, G = 22° [22| are closer to those 
in the FKS mixing scheme with = 43.7° and G = 33° 14J . 

In the rj-rj'-G mixing scheme, the situation becomes rather different for the decays 
B~ — > Gl~T>i with I = (e,/i, r), because of the effective mixing between the pseudoscalar 
glueball G and T] qiS . From the pQCD predictions for the branching ratios of Bu\ — > 
G(ll, lis, vis) as listed in Table EJ we observe that 

(i) The pQCD predictions for Br(B~ — > Gl~v{) with I = (e,/j,tau) are at the level 
of 10~ 6 , it is possible but not easy to observe such decay modes even at LHCb or 
future super-B experiments. 

(ii) The pQCD predictions for Br(B^ — > G(ll, lis, is is)) are sensitive to the mixing angle 
G , in fact they are proportional to sin 2 G according to Eq. lHBI) . as can be seen 
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TABLE IV. The pQCD predictions for the branching ratios of the considered decays in the 
rj-r]'-G mixing scheme. 



Decay modes 


<j> = 43.7° & 


' = 


33° 


<j) = 40° & G 




22° 


Br{B~ - 


-> rjl~vi) 


(0.331°;$ 


x 10" 


-4 


(0.37±8i*) 


x 


io- 


-4 


Br(B- - 


-)■ r}T~v T ) 




x 10" 


-4 


(o-23i8:8l) 


x 


10" 


-4 


Br(B~ - 


-> rfl-vi) 


(o.ie±85g) 


x 10" 


-4 


(o-i7l8:8S) 


x 


10" 


-4 


Br(B- - 


iT~V T ) 


(0.08i°;°I) 


x 10" 


-4 


(0.0918-g) 


x 


10" 


-4 


Br(B° - 


-> 7]l + l-) 


(o.39ig;ffl 


x 10" 


-8 


(0-45l8i 3 7 ) 


x 


10" 


-8 


Br(B° - 


-> r]T + T~) 


(0.80t°;l) 


x 10" 


-9 


(0-9218-1) 


x 


10" 


-9 


Br(B° - 


-» rjfu) 




x 10" 


-8 


(3-621^?) 


x 


10" 


-8 


Br(B° - 


■»■ r/i+r) 


(o-i9i8:8D 


x 10" 


-8 


CO 20 +0 - 08N | 


x 


10" 


-8 


Br(B° - 


-» T]'t + T~) 


(o.2o±g:8R 


x 10" 


-9 




x 


10" 


-9 


Br(B° - 


-> r/Vp) 


f 1 44+0-5i\ 

^•^-0.42^ 


x 10" 


-8 


(i-56±8:ti) 


x 


10" 


-8 


Br(B Q s - 






x 10" 


-7 


(2-2018:??) 


x 


10" 


-7 


Br{B° s - 






x 10" 


-7 


(o-48l8:? 8 ) 


x 


10" 


-7 


Br(B° s - 




/ 2 no+0.89^ 


x 10" 


-6 


/ 72 +0.76^ 


x 


10" 


-6 


Br(B Q s - 


-> r/i+r) 


(l-45±g:ig) 


x 10" 


-7 


^•^-0.67^ 


x 


10" 


-7 


Br(Bg - 


-» j/t+t - ) 


(0-18iHe 7 ) 


x 10" 


-7 


(o.24l8:S8) 


x 


10" 


-7 


Br(Bg - 


■> rjvv) 


(l-l<!o 7 ) 


x 10" 


-6 


(i.5ol8:l) 


x 


10" 


-6 



TABLE V. The pQCD predictions for the branching ratios of B Gl v\{t v t ) and — >• 
Gl + l~ (t + t~ , decays in the rj-r/'-G mixing scheme. 



Decay modes 


0G = 33° 


</>G = 22° 


Br(B~ -»■ Gl~T>{) 
Br(B~ GV - P T ) 


(0.64+8-1) x IO" 5 
(0.251*81) x 10-5 


(0.30l8iJ) x 10-5 
(0.12 ± 0.04) x IO" 5 


5r(i? -> Gi+r) 
Br(B° Gr+r-) 
J Br(5° Gi/i7) 


(0.7618:!') x io- 9 

(0.1818;8 5 7 ) x io -10 

(5.89+^t) x 10- 9 


(o.36l8;i?) x IO" 9 

(0.08+8-8^) x 10- 10 
(2.79lj;|i) x 10- 9 


5r(B s ° Gi+r) 
£r(£° Gr+r-) 
£r(.B° Gz^P) 


(0.2418;^) x lO-^ 
(0.88l8;|) x io- 9 
(1.8518:1) x io- 7 


(0.1ll8;8I) x io- 7 
(0.42l8il) x 10- 9 
(0.8818:1) x io- 7 



from the numerical results corresponding to different choice of 4>g — 33° or (pa = 22° 
in Table El 



(in) In Ref. [22| . the authors point out that 7^(1405) may be a possible candidate of the 
pseudoscalar glueball. We here define the ratio of the branching ratios Rg/ v and 
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find its value: 




BrjB- -+ Gl-vj) 
Br(B~ fjl-vi) 




0.19™, for0 G = 33°, 
0.081°$, for0 G = 22°. 



(47) 



If one can measure the ratio Br(B~ — > i](1405)l~vi)/Br(B~ — > r)l~i>i) experimen- 
tally, a direct comparison between the above pQCD prediction with the measured 
value will help us to determine whether the 77(1405) is the pseudoscalar glueball 
which have been searched for many years. 

V. SUMMARY 

In this paper we study the semi-leptonic decays B/B s — > (77,77' ,G)(l + l~ ,lv,vv) by 
employing the perturbative QCD factorization approach beyond the leading order. We 
firstly calculate the form factors of B/B s — > (r) q ,r) s ,r) g ) in the pQCD factorization ap- 
proach, show the g 2 -dependence of the form factors of B/B s — > (r) q ,r) s ,r) g ) transitions. 
And then we calculate the branching ratios of the corresponding semileptonic decay modes 
B/B s — >■ G)(l + l~ , lis, vv) in the conventional 77-77' mixing scheme and the r\-rj-G mix- 
ing scheme respectively. 

Based on the numerical calculations and the phenomenological analysis, We find the 
following points: 

(1) The NLO contribution to the form factors F O;+i t(0) for B/B s — > (rj q: r] s ) transitions 
can provide about 12% enhancement to the leading order ones. The NLO pQCD 
predictions for the relevant form factors i*o,+,r(0) agree well with the values obtained 
by using other methods. 

(2) For decays B~ — > r]l~9i and B~ — > r]'l~i>i with I = (e, //), the NLO pQCD predic- 
tions for their branching ratios agree well with the measured values. The pQCD 
predictions for Br(B~ — > rj^T~9 T ) are also at 10~ 5 level and could be measured in 
LHCb and forthcoming super-B experiments. 

(3) In the FKS mixing scheme, the pQCD predictions are Br(B Q s — > rj^vP) ~ 1.7x 10~ 6 , 
while in the range of 10~ 7 to 10~ 9 for other neutral decays. It is in fact a very hard 
job to observe these very rare neutral decay modes. 

(4) In the FKS mixing scheme, there is no mixing between the pseudoscalar glueball and 
the r] q and rj s . The branching ratios Br(B° — > Gl + l~) with I = (e, fi, r) are therefore 
smaller than 10 -11 , and can not be measured even at future super-B experiments. 

(5) The pQCD predictions for the branching ratios of the considered decay modes in the 
ordinary FKS mixing scheme are generally a little larger, about 20% in size, than 
those in the r)-r]'-G mixing scheme, except for the B® — > r]l + l~ (r + r~ , vv) decays. 

(6) In the rj-rj'-G mixing scheme, the pQCD predictions for the decay rates of the 
considered decays are sensitive to the value of 4>g- The pQCD predictions for 
Br(B~ — ?■ Gl~Vi) with I = (e, /i,r) are at the level of 10~ 5 , it is possible to observe 
such decay modes at LHCb or future super-B experiments. 
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(7) The pQCD prediction for the ratio Rc/rj is 



_ BriB -+GI vj) _ om 

Rc/r, - p /„_ 7— T - U.iy_ . 06 , (48) 

Br{B — > rjl vi) 



for (pc = 33°, which can be tested by the experiments. 
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Appendix A: Some relevant functions 

Here, we present the functions used in our calculation. The threshold resummation 
factors S t (x) is adopted from Ref. (23| : 

St(x) = 2l+2 ^l\ C) [x(l-x)r, (Al) 
V7rr(l + c) 

and we here set the parameter c = 0.4. The hard functions hi, h a , hb and h c come from 
the Fourier transform and can be written as 

h l (x 1 ,x 2 ,b 1 ,b 2 ) = KQ{^/x^X2~pmBbi)[0{bi - b2)Io( v ^x 2 ~pm B b 2 )K ( ^Jx2~pm B bi) 
+ 0(b 2 - bi)h{^Jx^m B b l )KQ(^fx2~pm B b 2 )}S t {x2)i 

h a (x 1 ,x 2 ,b 1 ,b 2 ) = K ( y /x 1 x 2 pm B b 2 )[9(b 1 - b 2 )I ( ^/xTpm B b 2 )K ( ^fx~[pm B b x ) 
+ Q(b 2 - b^Ioiy/xTpmsb^Koi y/xTpm B b 2 )}, 

h 2 (xx,x 2 , h, b 2 ) = 1 1 J-^_!_J^ B - 6{bx - b 2 )I (y/x^rjm B b 2 ) K (y/x^rjrrisbi) 

Z^/XiX 2 f]m B L 

+ Q{b 2 -b 1 )lQ(^fx^m B b 1 )K (^x^m B b 2 ) S t (x 2 ), (^2) 

h b (x 1 ,x 2 ,b 1 ) 



m 2 B xix 2 p - m 2 B [xi + (1 - xi)(l - x 2 )p] 
x [Ko{\/xi + (1 - £i)(l - x 2 )pm B bi) - K {y/xiX2pm B bi)], 
h c (x 1 ,x 2 ,b 1 ,b 2 ) = K (^xi + (1 -xi)(l -x 2 )pm B b 2 ) 

x [6>(&i - b 2 )h{\/x 1 (\ ~ P) + pm B b 2 )K (>y xi(l - p) + pm B bi) 
+ 0(b 2 - &i)/ (v^i(l - p) + pmsbijKa^xxil - p) + pm B b 2 )], 

where K , K\ and Iq are modified Bessel functions. 

The factors exp[— S B7] (t)] and exp[— S B c{t)] contain the Sudakov logarithmic correc- 
tions and the renormalization group evolution effects of both the wave functions and the 
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hard scattering amplitude with S Bv (t) = S B (t) + S v (t), where 

n / \ i m B j \ 5 f dp, . 

S B {t) = h) + o / — 7«(<MW), 

V2 d yi/bi 

m# , . -.rriB , N n [ da 

Sr,{t) = s{x 2 —7= ) b 2 ) + s((l - x 2 )^=, 6 2 ) + 2 / — 7 5 (a s (/x)) 
V2 V2 7i/6 2 A* 



(A3) 



(A4) 



with the quark anomalous dimension 7 g = — otJjK. The explicit expressions of the func- 
tions s(Q,b) can be found for example in Ref. [23]. The function Sg G (t) is defined as 



with 



and 



S BG (t) = S B (t) + S G {t), 



q u\ 1 P mB n, m \P mB h 
S G {t) = s G {x 2 —^, b 2 ) + s g {{1 - x 2 ) — — , b 2 j 



V2 



V2 



s G (Q,b) = 3 [ 
Ji. 



Q dp 
jb A* 



p 



71 



(A5) 
(A6) 

(A7) 



The hard scales t\ >2 and t a ^ c appeared in Eqs. (I25H30P are chosen as the largest scale 
of the virtuality of the internal particles in the hard 6-quark decay diagram, 

ti = max{ v /xipm jB , l/b 1} l/b 2 }, 

t 2 = t a = max{y/xTpm B , l/b x , l/b 2 }, 

t b = max{A/xi + (1 - Xi)(l - x 2 )pm B , ^x x x 2 pm B , I/61}, 

t c = max{ v / xi(l - p) + pm B , 1/bi, l/b 2 }. (A8) 



Appendix B: The effective Hamiltonian for b — > sl + l transition 
The effective Hamiltonian for the b — > sl~l + transition is of the form 



c 10 

* i=l 



where Ci(p) are the Wilson coefficients and are the local operators [37 

Oi = (s Q C a )y-A(c / 3& /3 )y-A, 2 = (s a Cp) V -A(c^b a ) V ^A, 

3 = (s a b a ) v „ A ^2(q'pq'p) v _ A , 4 = {s a bp) V - A ^%q' a )v-A, 



(Bl) 



5 = ( s Mv-A^2(q'pq' ) 
errib 



7 
9 



V+A 



O fi = (s n b 



07T Z 
^em 



8vr 



(j 7m j) [sy»(i - l5 )b) , o 



10 



aU0 V-A 



"87 



(W) [sy (1 - 75) b] 



(B2) 
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where q' = (u, d, c, s, b). 

amplitud* 

A(b sl+l 



The decay amplitude for b — > sl + l transition can be written as [37 

Gf a. 



-V*V tb {c l0 [s lfl (l- l5 )b] [lYlsl] 



2V2 vr 

+C 9 e// (/i) [ s - 7m (1- 75 )6] [IYI] 



-2m b C e 7 ff (fi) 



sia^y— (1 + 75) b 



where C 7 e// (/i) and C 9 e// (/i) are the effective Wilson coefficients, defined as 



Cf(fi) = CM + Y pcTt (s) + Y hB (s). 



(B3) 



(B4) 
(B5) 



Here the term Y p ert represents the short distance perturbative contributions and can be 
found in Ref. 38[. The term Yld(s), however, refers to the long-distance contributions 
from the resonant states and will be neglected because they could be excluded by experi- 
mental analysis [3 9l. |4Cfl. The term C' b ^ in Eq. ( 1B4|) is the absorptive part of b — > 37 and 
is given by [4 If 



-r/ 14/23 [Giixt) - 0.1687] - 0.03C 2 (/x) 



with 



GAx t ) 



%t \ x t 0J/ t 



2) 3x 2 t \nx t 



where r\ = a s {m w ) j a s (fi) and x t = m^jrn^. 



+ 



4(x t -l) 4 ' 



(B6) 



(B7) 
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